Effect of a columnar defect on the shape of slow-combustion fronts 
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We report experimental results for the behavior of slow-combustion fronts in the presence of a 
columnar defect with excess or reduced driving, and compare them with those of mean-field theory. 
We also compare them with simulation results for an analogous problem of driven flow of particles 
with hard-core repulsion (ASEP) and a single defect bond with a different hopping probability. The 
difference in the shape of the front profiles for excess vs. reduced driving in the defect, clearly 
demonstrates the existence of a KPZ-type of nonlinear term in the effective evolution equation for 
the slow-combustion fronts. We also find that slow-combustion fronts display a faceted form for 
large enough excess driving, and that there is a corresponding increase then in the average front 
speed. This increase in the average front speed disappears at a non-zero excess driving in agreement 
with the simulated behavior of the ASEP model. 
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I. INTRODUCTION 

Nonequilibrium interfaces that display interesting scal- 
ing properties are quite common in physical (crystal 
growth, fluid penetration into porous media, etc.), chem- 
ical (reaction fronts), as well as biological (growing bac- 
terial colonies) systems. The dynamics of these systems 
have long been thought to be generically described by the 
Kardar-Parisi-Zhang (KPZ) equation jjj, or some other 
equation of motion in the same universality class 0] . In 
two space dimensions in particular (one-dimensional in- 
terfaces) when exact solutions are available, the scaling 
properties of the KPZ equation are well understood. The 
same is not, however, true for the experimental observa- 
tions of scaling of interfaces. Typically one has found 
a roughening exponent clearly higher than that for the 
KPZ equation 0. Various explanations have been sug- 
gested as for why the KPZ scaling has not generally 
been found, and most of the time correlated and/or non- 
Gaussian noise have been the prime suspects 0, ■ 

Recent experiments on slow-combustion fronts propa- 
gating in paper 0, || , and on flux fronts penetrating a 
high-T c thin-film superconductor have provided new 
insight into this problem 0. It indeed appears that 
short-range-correlated noise, quenched and dynamical, 
with possibly at the same time a non-Gaussian ampli- 
tude distribution for small time differences, induce an 
additional time and an additional length scale, beyond 
which KPZ scaling can only be observed. Despite these 
recent advances, it would still be worth while to demon- 
strate the existence of the nonlinear term, as introduced 
by Kardar, Parisi and Zhang ||, and essential for the 
KPZ dynamics, directly from the observed fronts. This 
would in essence prove that KPZ type of dynamics, in- 
cluding possibly effects of nontrivial noise, can indeed 



be expected to be generic for nonequilibrium interfaces. 
One can demonstrate the presence of this term indirectly 
by using, e.g., inverse schemes able to infer the (partial) 
differential equation that governs the observed stochas- 
tic evolution of interfaces [IJ , but there is also a way to 
produce a directly observable effect on the shape of the 
interface, due to this term. 

This method for observing the operation of the non- 
linear term was suggested already for some while ago by 
Wolf and Tang . They considered the effect of colum- 
nar defects, columnar in two space dimensions on which 
case we concentrate here, and found that there is a clear 
'asymmetry' between the shapes of the fronts that prop- 
agate in the presence of an advancing (excess driving) 
and a retarding (reduced driving) defect. This asymme- 
try is a direct consequence of the nonlinear term in the 
KPZ equation. For a positive coefficient in this term, 
applicable to slow-combustion fronts, the noise-averaged 
front should be faceted with a forward-pointing trian- 
gular shape around an advancing defect, with a height 
proportional asymptotically to the width of the sample 
in the case of a defect in the middle of the sample, or to 
the basic period in the case of periodic boundary condi- 
tions. In the case of a retarding defect, the shape of the 
front should not be faceted, and the (negative) height of 
the deformation in the profile should be proportional, ac- 
cording to this mean-field theory, to the logarithm of the 
basic period. Despite its apparent simplicity, this kind of 
experiment has never been performed. 

As is well known, the BCSOS interface model in which 
the nearest-neighbor heights are restricted to differ only 
by ±1, displays KPZ behaviour, and is on the other 
hand equivalent to a driven flow of particles (hopping 
ratep) with hard-core repulsive interactions (ASEP) @|. 
A columnar defect in an interface model corresponds to 
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a fixed slow or fast bond (hopping rate rp with r < 1 
or r > 1, respectively) in the ASEP model. A faceted 
interface corresponds to a traffic jam of infinite length 
in the thermodynamic limit behind the slow bond. An- 
other related question is the detailed shape of the den- 
sity/interface profile. The mean-field theory of ^| pre- 
dicts an infinite queue for all r < 1 and no queue (a 
logarithmic decay of the density profile) for r > 1, i.e., 
r c = 1. Janowsky and Lebowitz 9] considered the to- 
tally asymmetric TASEP model with a slow bond, but 
concentrated mainly on the shock wave fluctuations far 
away from the slow bond, and apparently did not con- 
sider the faceting/queueing transition (they had a phase 
diagram with r c = 1). The same model was also con- 
sidered by Kolomeisky [lQj, but he did not consider the 
faceting/queueing transition either. 

Kandel and Mukamel considered a somewhat different 
model, which is sup pos ed to be in the same universality 
class, and proposed |14j that the faceting/queueing tran- 
sition should take place at an r c < 1. Their simulation 
data were not, however, conclusive. 

In slow-combustion experiments the detailed shape of 
the front profile is difficult to determine, and thereby also 
the disappearance of faceting. Faceting is however re- 
lated to increased front speed, also in the thermodynamic 
limit, and this is an easier observable. For possible non- 
faceted fronts, which would correspond to r c < r < 1, an 
increased front speed would only be a finite-size effect, 
as also the decreased front speed in the case of a retard- 
ing defect corresponding to r > 1. Notice that the ef- 
fective nonlinear term is positive in the slow-combustion 
experiments, while it is negative in the ASEP models. 
Therefore, an advancing (retarding) columnar defect in 
the first case corresponds to a slow (fast) bond in the 
latter case. 



II. EXPERIMENTAL DETAILS 

The equipment we use in slow-combustion experiments 
has been described elsewhere 0, , so it suffices to say 
here that samples were 'burned' in a chamber with con- 
trollable conditions and that the video signal of propagat- 
ing fronts was compressed and stored on-line on a com- 
puter. The spatial resolution of the set up was 120 /im, 
and the time resolution was 0.1 s. For the samples we 
used the lens-paper grade (Whatman) we have used also 
previously [j| . Lens paper was now used to speed up the 
experiments. 

As slow-combustion fronts do not propagate in paper 
without adding an oxygen source for maintaining the 
chemical reaction involved, we added as before [j, |5| a 
small amount of potassium nitrate in the samples. This 
method also allows for a relatively easy way to produce 
advancing and retarding columnar defects. By using 
masks it is straightforward to produce a narrow (ver- 
tical) stripe with a smaller or an additional amount of 
potassium nitrate. The average concentration of potas- 



sium nitrate determines the average speed of the fronts, 
so it serves as the control parameter of the problem. 

It is, however, quite difficult to accurately regulate 
the amount of potassium nitrate absorbed in the sam- 
ple. This means that it is difficult to produce samples 
with exactly the same base concentration, and the same 
concentration difference between the base paper and the 
columnar defect. Therefore, the statistics we get for any 
fixed difference in the concentration is not quite as good 
as we would hope. They are adequate for the main fea- 
tures of the fronts but not for such details as, e.g., ac- 
curate forms of the front profiles. They are also good 
enough for a quantitative analysis of changes in the front 
speed. 

The samples were typically 20 cm (width) by 40 cm, 
and the columnar defect (vertical stripe) in the middle 
of the sample was 1.0 cm wide. The defect cannot in our 
case be too narrow as fluctuations in the slow-combustion 
process would then tend to wipe out its effect. Too wide 
a stripe would on the other hand cause effects due to its 
nonzero width, which are unwarranted. We also used 
simulations with a discretized KPZ equation to check 
that the ratio 1 cm to 20 cm should not cause addi- 
tional effects [l^. The length of the samples was in 
most cases adequate for achieving stationary behavior, 
and only those results are used here where saturation of 
the profile was evident. 

When analyzing the front profiles, the stripe was re- 
moved from the data, as well as about 6 mm from both 
boundaries of the samples. As the system is symmetric 
across the stripe in the middle, the observed front profiles 
were also symmetrized for better statistics. 

As reported already before 0, fluctuations in the 
slow-combustion fronts in paper are noticeable. For ex- 
treme values of the potassium-nitrate concentration there 
appear problems with pinning (low concentration) or lo- 
cal avalanche type of bursts (high concentration) in the 
fronts. Also, very small values of the concentration differ- 
ence between the base paper and the defect stripe could 
not be used, as fluctuations then completely masked the 
effect of the defect. These problems were noticeable for 
retarding defects in particular. In the data reported here, 
concentration varied between 0.265 and 0.61 gm" 2 in the 
base paper, between 0.1 and 1.05 gm -2 in the stripe, and 
the concentration difference varied between 0.06 and 0.49 
gm on the positive side (25 burns), and —0.197 and 
—0.478 gm -2 on the negative side (19 burns). Because 
of the practical restrictions and fluctuation effects, the 
number of successful burns was relatively small. 



III. RESULTS 

We will need the dependence on potassium-nitrate con- 
centration of the front velocity below so we consider it 
first. It is useful to begin with a discussion of the accu- 
racy of the front-velocity determination. 
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FIG. 1: Velocity distributions for slow-combustion fronts with 
potassium-nitrate concentrations 0.34 gm~ 2 (full line), and 
0.536 gm -2 (dashed line). Velocities are determined for a 
time difference of 2 s. 



FIG. 2: Front velocity as a function of potassium-nitrate con- 
centration. Measured points are denoted by circles and the 
line is a linear fit to these points. 



A. Front velocity 

Lens paper is thin so that variations in its mass density 
and dynamical effects such as (possibly turbulent) con- 
vection around the combustion front are both expected 
to give a contribution to the effective noise. Noise ampli- 
tude is consequently relatively large [4, 5], and therefore 
also velocity distribution of a propagating front can be 
expected to be broad. We show in Fig. 1 that distri- 
bution for two different values of the potassium-nitrate 
concentration. 

It is evident from this figure that, even though the 
average velocities in these two cases are relatively well 
separated and easily distinguishable, the velocity distri- 
butions have a big overlap. Together with the limited 
statistics for any fixed value for (the difference in) the 
potassium-nitrate concentration, these broad distribu- 
tions mean that some variation can be expected to occur 
in the measured average velocities (velocity differences) 
of the fronts. 

In the presence of a columnar defect, we should, e.g., 
determine the change in the average front speed arising 
from the defect. This can be accomplished by analyzing 
separately the undeformed 'flat' part of the fronts, and 
the part of the front profile affected by the presence of 
the defect. Determination of the average speed of the fiat 
fronts is done in the transient (wrt profile shape) phase 
in which the (growing) width of the deformed profile is 
still less than the width of the sample. In this phase the 
fiat part of the front is already in the saturated regime 
with constant average velocity. The average speed of the 
deformed profile is determined in a later phase in which 
the width of the deformed part of the profile essentially 
coincides with the sample width. 

We have also determined the average front velocity 
for 122 individual burns for a fairly broad interval in 



the potassium-nitrate concentration, and these data are 
shown in Fig. 2 together with a linear fit to the measured 
points. 

The dependence on potassium-nitrate concentration of 
the front velocity is not expected to be linear especially 
near the pinning limit, but, for the concentration range 
shown here, it is well approximated by a linear behavior, 
which is also more convenient for the subsequent analy- 
sis. We find that, with a linear fit to the data, the front 
velocity v is given on the average by 



4.2C + 6.2, 



(1) 



where C is the potassium-nitrate concentration, and v is 
in units mms -1 when C is expressed in gm . 

Before showing the measured front profiles in the pres- 
ence of a columnar defect, let us consider, in order to 
make later comparisons more transparent, what is ex- 
pected from the mean-field solution as reported in Ref. 



B. Mean-field prediction 

We assume that the time evolution of the fronts h(x,t) 
is governed by the KPZ equation dh/dt — vV 2 h + 
■|(V/i) 2 + k + r)(x, t), where r](x, t) describes white noise 
with delta-function correlations in space and time, and 
the driving term contains the idealized defect as a delta- 
function contribution, k — kq + k± y\ 5 (x — L/2 + nL). 
For simplicity we assume here as in [hit periodic bound- 
ary conditions. If we average over noise in the KPZ equa- 
tion, and denote H(x,t) = (h(x,t)), we find that 



dH(x,t) 

Ft 



vH" + i(H>) 2 +'i((V5hf)+K, (2) 
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in which H' = H'(x,t) denotes the spatial derivative 
of H, and Sh = h(x, t) — H(x, t) describes fluctuations 
around the noise-averaged profile. A corresponding equa- 
tion can be derived for Sh |l3[ . 

As Sh should not depend (locally) on H(x,t) nor 
on H'(x,t), and only local interdependence between Sh 
and the noise averaged profile can be assumed to ap- 
pear, one would then expect |13| that in leading order 
(CVSh) 2 ) — clq + a,2H"(x,t), with ao and 02 some con- 
stants. This assumption will make Eq. J2J) closed so that 
it can be solved without further reference to the fluctu- 
ations. The delta- function contribution in the driving 
term will induce cusps in H(x) at x = L/2 — nL 1 and the 
solution of Eq. (0 is therefore equivalent to solving the 
equation 

d -^ = v e H" + \{H'f + n e (3) 

in the interval —L/2 < x < L/2, with boundary condi- 
tions H(L/2) = H{-L/2) and H'(±L/2) = ±s. Here 
we have defined the effective (renormalized) parameters 
v e = v ± Xa 2 /2 and n e = k q + Aao/2, and the magnitude 
of the slope of the front at the defects is s = Ki/2v e . 

Equation © is the well-known Burgers equation 
which can be solved in closed form in one space dimen- 
sions. It is useful to express it first in dimensionless 
form, which can be achieved with transformations H = 
H H, x = H Q x/s, t = 2H t/(\s 2 ), with H = 2v e /\ the 
internal length scale of the system. 

We look for a stationary solution of this equation in 
the form H{x,t) = {n e + sgn(Ki)q 2 )t ± ln(f(qx)), where 
sgn(,z) is the sign of z, and we have already used the 
Hopf transformation in the spatial part of the ansatz to 
remove the nonlinearity from the equation for /. We find 
that 13] 



f(z) = cosh(z), gtanh ^— j = 1, (4) 
for advancing defects (ki > 0), and 



f(z) = cos(z), gtan ^— j = 1, (5) 

for retarding defects (ki < 0). Here L = sL/Hq is 
the dimensionless width of the system, and z = qx. 
Asymptotically, for L ^ 1 (and A > 0), the profile 
around an advancing defect is a forward-pointing tri- 
angle with sides that have slopes ±s, and with height 
AH + = H(L/2) - H(0) ~ si/2. The asymptotic profile 
around a retarding defect is given by -ffoln|cos(7ra:/Z/)| so 
that Ai7_ ~ ~H \n(sL/nH ). The magnitude of H + 
thus grows linearly with L (or s) while that of only 
grows logarithmically with L (or s). This asymmetry is 
a direct consequence of the nonlinear term that enhances 
the deformation in the former case but reduces it in the 
latter case. 
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FIG. 3: Average front profiles with a columnar defect for 
AC = ±0.33 (upper and lower profile, respectively). 



C. Measured front profiles 

In the above mean-field theory, 8v = k\/L is the dif- 
ference between the front velocity (driving) inside the 
defect and outside the defect. In the slow-combustion 
experiments this velocity difference is regulated by 
the potassium-nitrate concentration so that now Sv — 
4.2AC, where the numerical factor comes from the linear 
fit given by Eq. Q , and AC = Cdefect — Cbasc is the con- 
centration difference. This means that the scaling factor 
s is given by s = 4.2L\AC\/v e . Without as yet knowing 
the actual value of v e needed for evaluating the size of s 
and Hq, reasonable estimates, based on the results from 
the inverse method solution for the effective equation of 
motion [12| , indicate that the slow-combustion fronts are 
not necessarily in the strictly asymptotic regime: we ex- 
pect that L > 1 but not by a very big margin. Notice 
that the size of L is now regulated by AC as the width L 
of the samples is held fixed. Despite the achievable val- 
ues of L, we can expect to clearly see the 'asymmetry' in 
the heights of the front profiles for different signs of the 
concentration difference. In Fig. 3 we show the averaged 
(and symmetrized) front profiles for AC = ±0.33. 

It is indeed evident that there is a clear difference in 
the heights of the front profiles around advancing and re- 
tarding defects. By following the transient time evolution 
of the fronts, we could also see a clear difference there. 
For AC > 0, when a triangular deformation was formed 
after a while around the central stripe, its height and base 
length grew with a more or less constant velocity until 
the base length reached the width of the sample, while 
the slopes of the sides of this triangle remained roughly 
constant. For AC < on the other hand, the height 
of the deformation saturated much faster even though it 
also grew more or less linearly in time in the beginning, 
and the base length of the deformation reached the sam- 



5 



600 




FIG. 4: Successive fronts with a time difference of 0.5 s for 
concentration difference AC = 0.327. Also marked are the 
stripe with the excess concentration of potassium nitrate, the 
fronts between which the average profile is determined (thick 
lines), the height of the final profile, and the average shape of 
the profile. 



pie width at the same time. This transient behavior will 
be analyzed in more detail below. The qualitative behav- 
ior for the AC > case is clearly visible in Fig. 4 which 
shows the successive fronts with a time difference of 0.5 
s for AC* = 0.327. 

A more quantitative comparison between the mean- 
field solution for the noise averaged front and the ob- 
served slow-combustion fronts can also be made. For 
this purpose we have found it convenient to consider in- 
stead of the profile heights AH± the average slopes of 
the left-hand (LH) sides of the profiles (c.f. Fig. 3), 
k± = 2AH±/L. As we do not expect to be in the strictly 
asymptotic regime, we have used the full transcendental 
equations for q in Eqs. and JSJ) above when fitting 
the observed k± with the mean-field result. 

The average slopes, as functions of concentration dif- 
ference AC, will now depend on two parameters, A = 
Hq/L and B = 2.1L/v e , which are used to fit the mea- 
sured slopes. From the fitted values for these parameters 
we can then estimate the coefficients A and v e for this 
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FIG. 5: The absolute value of the average LH slope of the 
deformed front profile around a columnar defect as a function 
of difference in the potassium-nitrate concentration. The full 
line is a fit by the mean-field solution with Eqs. and Q). 



system. 

We show in Fig. 5 the experimentally determined val- 
ues for k + and fc_ together with the fit by the mean-field 
solution using Eqs. @ and JSJ. 

Fits to the data were not very sensitive to the actual 
value of parameter A so that the correlation coefficient 
did not change much even if A was changed in a relatively 
large interval. If the AC > and AC < data were fit- 
ted separately without any restrictions on the two param- 
eters, these fits had also a tendency to produce somewhat 
different values for the two cases. As the signal-to-noise 
ratio is better for the AC > data, we fixed A such that 
it was between the two separately fitted values but closer 
to the one from the unrestricted two-parameter fit to the 
AC > data, and in the interval within which the quality 
of this fit was essentially unchanged: A ~ 0.3. Thereafter 
an unrestricted one-parameter fit to the whole data was 
used to find the value for B. In this way we found that 
B ~ 2.5. 

The fitted values for parameters A and B allow now an 
estimation of two physical parameters, the 'renormalized' 
diffusion coefficient v e and the coefficient of the nonlinear 
term, A. We thus find that v e ~ 144 rrrm. 2 s _1 , and A ~ 
5.6 mms -1 . In the estimate for v e we used an 'effective' 
sample width L e g ~ 180 mm, a bit smaller than the 
width 202 mm of the actual sample, due to the width 
of the defect stripe and to allowing for some boundary 
effects. By other methods we have found previously that 
A ~ 4.1 — 5.1 mms -1 [l^, so that the value found here is 
fairly close to these previous estimates. 

In view of the unavoidable fluctuations in the measured 
averaged slopes, we find the fits to the measured points 
by the mean-field solution to be quite reasonable. 
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D. Defect-induced change in front velocity and the 
queueing transition 

As already discussed above, the mean-field solution 
predicts a faceting or queueing transition at AC = 0. 
Above this transition (AC > 0), the average front ve- 
locity is increased due to the presence of an advancing 
columnar defect, and below this transition the change in 
front velocity should vanish for large enough L. For neg- 
ative AC the change in velocity is negative, and should 
decrease in magnitude with increasing |AC|. According 
to Kandel and Mukamel 0], this transition should ap- 
pear at a AC cr > 0. 

As the numerical data of [1J| is not decisive, we have 
done simulations on a totally asymmetric ASEP 
model with a fixed defect bond with hopping rate rp in 
the middle of the system, while the hopping rate at the 
other bonds was p. Open boundary conditions were im- 
posed such that the hopping-in rate at the left boundary 
was ap, and the hopping-out rate at the right bound- 
ary was Pp. In what follows we only consider the case 
a = 0=p=l/2. 

This model shows |ll| a queueing transition at r = 
r c = 0.80 ± 0.02. In addition, the density profile dis- 
plays a qualitatively similar asymmetry between the slow 
and fast defect-bond cases as the mean-field solution 
for the KPZ fronts between the advancing and retard- 
ing columnar-defect cases. In the thermodynamic limit 
the deformation stays non-zero only in the faceted phase 
above the transition. The density profiles in both the 
faceted and non-faceted phases also display interesting 
power-law tails. 

As the detailed shapes of the front profiles are difficult 
to determine experimentally, we only compare the results 
for the dependence of the average front velocity V = (v) 
(current J = (j)) on the potassium nitrate concentra- 
tion C (hopping rate p). In this comparison dimension- 
less variables are used, (V — Vq)/Vq for the change in 
the front velocity, and similarly for the current but for 
reversed sign as an advancing defect corresponds to a 
slow bond, AC /Co for the potassium-nitrate concentra- 
tion difference, and Ap/p = 1—r for the hopping-rate dif- 
ference. Differences are all determined between the value 
with or at the defect and the value elsewhere or without 
the defect. In this way no fitting is involved in the com- 
parison. Obviously the actual driving force is not known 
exactly for the slow-combustion fronts, but the observed 
linear dependence well above the pinning transition be- 
tween the potassium-nitrate concentration and the front 
velocity suggests that the dimensionless difference can be 
reliably used in this kind of comparison. 

This comparison of the slow-combustion experiment 
and the totally asymmetric ASEP model results is shown 
in Fig. 6. It is evident from this figure that agreement 
between the two results is reasonable as there is no fit- 
ting involved. There are still fairly large fluctuations in 
the experimental data, and it is not possible to have re- 
sults for very small values of AC as fluctuations tend to 




(C-C )/C , 1-r 

FIG. 6: Scaled velocity change of slow-combustion fronts due 
to a columnar defect as a function of (C — Co) /Co (*), and 
scaled current change in the totally asymmetric ASEP model 
due to a defect bond as a function of (p — rp)/p = 1 — r (o). 
The full line connects the latter data points. 

wipe out the whole effect, and the system is then not in 
the 'asymptotic regime'. These results indicate, however, 
that there indeed is a faceting (queueing) transition at a 
non-zero AC cr (Ap cr , i.e., r c ^ 1). 

E. Transient behavior 

In addition to the stationary profiles analyzed above, 
it is also possible, as already indicated, to study the tran- 
sient profiles, i.e., how the defect-induced profiles grow 
at the initial phases of the process. The transient be- 
havior of the profile around an advancing column is par- 
ticularly simple. The Burgers equation Eq. Q admits 
in this case a solution of exactly the same shape as the 
stationary solution, which grows linearly in time until its 
baseline reaches the width of the sample. Such a 'self- 
similar' transient does not exist in the case of negative 
AC, so analytical results for transient behavior are then 
difficult to find. In the non-asymptotic regime L « 1 
one can however show that the situation is symmetric, 
A_ff_ ~ — AH + . One would thus expect that, at least 
in our case when L is not particularly large, the height 
|AiJ_| would also grow initially (at least nearly) linearly 
in time. 

The expected transient behavior for AC > is already 
(qualitatively) evident from Fig. 4 above. More quanti- 
tatively the transient time evolution of the height of the 
deformed profile can be analyzed, e.g., by plotting H(2t) 
against H(t). For a linear time evolution the former value 
is twice the latter. In Fig. 7 we show this plot, averaged 
over 32 individual burns, including both signs of AC. 

The initial transient behavior is approximately linear 
in time for both cases. For AC > the trend continues 
nearly linear until saturation sets in when the width of 
the profile equals the width of the sample. For AC < 
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H(t) [mm] 

FIG. 7: AH(2t) as a function of AH(t) averaged over 32 
burns. Positive values correspond to AC > and negative 
values to AC < 0. The full line is AH(2t) = 2AH(t). 

the behavior is quite similar except that saturation takes 
place earlier. There is also some indication that, in this 
case, the growth of AH becomes nonlinear in time al- 
ready before saturation, but the quality of the data does 
not allow for a decisive conclusion on this. 



IV. DISCUSSION AND CONCLUSIONS 

The difference in the amplitude (height), and perhaps 
not so clearly in the shape, of the front in the slow- 
combustion experiments, caused by a columnar defect 
with excess or reduced driving, respectively, was clearly 
demonstrated. The behavior of the height of the de- 
formed profile, and the qualitative shape of the profile 
in the case of excess driving, were also reasonably well 
explained by the mean- field solution of Ref. fTs]. The 
asymptotic shape of the profile in the case of negative ve- 
locity difference could not be unequivocally determined 
as fluctuations are more important in this case of rela- 
tively small amplitude of the profile. The reduced, in 
comparison with the case of excess driving, height of the 
profile was very evident. In the case of positive velocity 
difference the transient behavior of the profile, i.e., the 
growth of the defect induced deformation in the profile 
shape, could as well be explained by the mean-field so- 



lution. For negative velocity difference a nearly linear 
behavior in time was observed initially, followed perhaps 
by a regime of nonlinear time evolution before saturation. 

Fitting the average height (or equivalently the average 
slopes of the sides) of the profile with the mean-field so- 
lutions provided us with estimates for the effective 'diffu- 
sion constant' v e and the coefficient of the nonlinear term, 
A. The latter parameter can also be determined from the 
slope dependence of the local front velocity 0, 0] , or by 
applying an inverse method on the observed fronts |12| . 
The value found here for A is fairly close to these previ- 
ous estimates, and we find this level of agreement very 
reasonable in view of the rather large fluctuations in the 
present data. 

One should, however, notice that the A measured for a 
sample depends on the potassium-nitrate concentration 
in that sample, and that the average potassium-nitrate 
concentration was not the same in the samples used in the 
experiments. We did not take this variation into account, 
as it can be assumed to give a small effect in comparison 
with the other experimental uncertainties, so that the 
present estimate represents an 'average' value. 

The effective diffusion coefficient v e contains, in ad- 
dition to the bare diffusion coefficient of the original 
KPZ equation, an unknown rcnormalization factor due 
to noise-induced fluctuations around the average front 
profile. We cannot thus get an estimate for the 'bare' 
diffusion coefficient which can be estimated by other 
means [l^. However, we can conclude that the noise- 
induced renormalization of v & appears to be sizable. 

The position and nature of the faceting (queueing) 
transition in interfaces affected by a columnar defect (in 
the ASEP model by a defected bond), has been a long- 
standing problem. The agreement found here between 
slow-combustion experiments with a columnar defect and 
the related TASEP model results, indicates that this 
transition is indeed at a non-zero value of the respective 
control parameter. No scaling properties of the transition 
could be analyzed at this stage, but the TASEP model 
results also indicate that this transition is continuous. It 
remains an experimental challenge to analyze this tran- 
sition in more detail. 
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